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The differential rates and emissivities of neutrino pairs from an equilibrium plasma are calculated 
for the wide range of density and temperature encountered in astrophysical systems. New analytical 
expressions are derived for the differential emissivities which yield total emissivities in full agreement 
with those previously calculated. The photon and plasmon pair production and absorption kernels 
in the source term of the Boltzmann equation for neutrino transport are provided. The appropriate 
C**) • Legendre coefficients of these kernels, in forms suitable for multi-group flux-limited diffusion schemes 

' are also computed. 
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I. INTRODUCTION 



A wide range of astrophysical phenomena are significantly influenced by weak interaction processes that involve 
the emission or absorption of neutrinos in matter at high density and/or temperature. Examples include neutrino 
energy loss in degenerate helium cores of red giant stars |2| , cooling in pre- white dwarf interiors [3j , the short- and 
long-term cooling of neutron stars 0, Q , the deflagration stages of white dwarfs which may lead to type la supernovae 
' 0, 13 1 explosive stages of type II (core-collapse) supernovae Q , and thermal emission in accretion disks of gamma-ray 
bursters pi llflj. (The selected references contain more complete references to prior and ongoing work.) Depending 
on the density and temperature of ambient matter, the emission of neutrinos via both neutral and charged current 
interactions is an important energy-loss mechanism, while scattering and absorption of neutrinos serve to deposit 
energy into matter. In recent years, it has been realized that neutrino cross-sections in matter and the implementation 
of accurate neutrino transport are critical to understand the precise mechanisms which trigger explosive events. A 
sterling example is provided by gravitational-core-collapsc supernova simulations in which a strong coupling between 
neutrino transport and hydrodynamics (in many cases supplemented by convection, rotation, and magnetic fields) 
| has been realized [Tlj . 

Neutrino transport in the supernova environment is described by a Boltzmann transport equation, which is a 
nonlinear integro-partial differential equation that describes the time rate of change of the neutrino distribution 
function /. Advances made to date in the numerical solution of this equation in the supernova context may be 
found in Refs. Historically, multigroup methods (in which the equation is discretized in energy groups) have 

involved the use of moment equations. When the temporal derivative of the first order deviation in / is set to zero, 
a diffusion equation is obtained, but this cannot adequately handle the free-streaming regime at low densities. Flux 
limiting schemes have been used to bridge the diffusive and free-streaming regimes, but these are somewhat arbitrary 
and calibrations vary depending upon neutrino opacities and dynamics. In addition, there is a coupling of different 
neutrino-energy groups, especially because of neutrino-electron scatterings which involve large energy transfers. An 
additional complication in supernovae is that the approach to thermal and chemical equilibrium, and the passage 
from diffusive flow to free streaming, require careful treatment. Even with modern parallel supercomputers, it is 
necessary to integrate the Boltzmann equation over solid angles to reduce the dimensionality of the problem, with 
a corresponding loss of information about the neutrino angular distribution function. This could be important in 
regimes in which neutrino-driven convection, a 3-D phenomenon, is occuring. 

The basic microphysical inputs of accurate neutrino transport coupled in hydrodynamical situations are the dif- 
ferential production and absorption rates and their associated emissivities. The precise forms in which such inputs 
are required for multienergy treatment of neutrinos is detailed in Ref. ^2 ( see a l so R- e f- US m which other current 
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developments are summarized). The objective of this work is to make available differential rates and emissivities 
for the thermal production of neutrino pairs from a plasma (7* — > v + v) . Similar quantities for the photoneutrino 
process, e + ^^e + v + D, will be presented separately. 

We wish to note that the total neutrino pair emissivity from the plasma process has been investigated in several 
prior works [3 El El > an d a comprehensive treatment of this process valid over a wide range of density 

and temperature can be found in Ref. [T^j . (The results of this latter reference are used extensively in this work.) 
However, in prior works in which total emissivities were computed the energy and angular dependences of the emitted 
neutrinos were eliminated with the help of Lenard's identity: 

/ ^^ 4 (K-Q 1 -Q 2 )Q^ 2 = ^e(K 2 )(2K^K^ + K 2 g n, (1) 

where = (E\,qx) and Q2 = (£ , 2,q2) are the 4-momenta of the outgoing neutrinos, and A' M = (u;,k) is the 
4-momentum of the massive photon in the medium. While the use of Eq. JQ) simplifies considerably the calculation of 
the total emissivity, differential information about the neutrinos is entirely lost. In addition, calculations of differential 
rates and emissivities entail the calculation of the relevant squared matrix elements, which was hitherto bypassed in 
obtaining the total rates and emissivities. 

Section [H] is devoted to obtaining the squared matrix elements which are employed in the calculation of the 
differential rates and emissivities discussed in detail in Sec. 11111 A check of these results is provided by comparing 
the total emissivities obtained by both analytical and numerical integrations of the differential emissivities. The 
results for the total emissivities reported in Ref. [3] are used as an additional check of our results. Qualitative and 
quantitative discussions of the numerical results obtained here are presented in Sec. II VI Sec. [V] contains production 
and absorption kernels in forms suitable for detailed calculations of neutrino transport. In order to assess the relative 
importance of neutrino pair production from the plasma process, a comparison with competing processes is made 
in Sec. IVII over a wide range of density and temperature. A summary of this work is provided in Sec. IVIII For 
completeness, a brief description of the standard model effective coupling is included in Appendix [S] Except when 
presenting numerical results, we use units in which h = c = ks = 1- 
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FIG. 1: Leading order Feynman diagrams describing the emission of a neutrino pair from the plasma process. The charged 
current process in which the VF-boson is exchanged produces only v e v e , while that in which the neutral Z-boson is exchanged 
results in pairs of all three neutrino (e, fx, and r) flavors. 



II. PHOTON AND PLASMON DECAY 



Our discussion will be restricted to the case of an equilibrium plasma in which the net negative electric charge 
of electrons and positrons is cancelled by a uniform positively charged background of protons, alpha particles, and 
heavier ions. The equation of state and the phase structure of matter, and the abundances of the various constituents 
including those of dripped neutrons at sub-nuclear densities are determined by the minimization of free energy. 

As is well known, e + e~ pairs in a plasma cause the photon to acquire an effective mass K 2 = II ^ 0, which arises 
from electromagnetic interactions (cf. Ref. [T^j and references therein). Therefore, we can consider the photon to 
be a massive spin-1 particle that couples to the vv pair through the two one-loop diagrams shown in Fig. ^ The 
channel in which the exchange of the Z-boson occurs can produce any of the three species (ve.^.r) of neutrinos and 
their anti-particles, whereas the channel in which the IF-boson is exchanged, only the v e v e pair is produced. 
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A. Effective Coupling 

At the low energies of interest here, the matrix element for neutrino pair production from the plasma is given by 



u 1 {Q l )Y{l-lMQ2) 



(2) 



where is the photon's polarization, T^ V {K) is the effective vertex tensor, u\ and V2 are the wave functions of the 
outgoing neutrinos, and the Dirac 7-matrices have their usual meaning. The 4-momcnta A", Qi and Q2 arc defined 
earlier in connection with Eq. The effective vertex can be expressed as 
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where e and m e are the charge and mass of the electron, Gf is Fermi's weak interaction coupling constant, and T is 
the temperature. C v and G A are the vector and axial couplings for a neutrino of flavor /. The trace is taken over 
the spin states and denotes the sum over the Matsubara frequencies 



P = (2/ + l)7iTi + [i e 



(4) 



of the e + e loop in a heat bath with chemical potential [i e . The e and e + loop momenta are denoted by Pf = 
(E el ,pi) = P^ = (E e ,p) and Pg = (E e2 ,p 2 ) = P' 1 - K^, respectively. The T^K) in Eq. © satisfies current 
conservation, i.e., K^T^iK) = 0. It is advantageous to decompose the vertex tensor into its vector and axial 
parts: 



T. V (K) = V. v {K)+iA. v {K), 
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(5) 



Summing over the Matsubara frequencies and using contour integration, V^ V (K) and A^ V (K) can be expressed as 



ir / tj-\ , eG F C v f d 3 p {K»P V + KyPjjK ■ P) - K 2 P^P U - g^(K ■ P) 2 / r^,- (w 7^ 

V(A) = 4 72W J (A ■ - (A 2 ) 2 /4 (n F (P e ,M e ,r) + n F (^,Me,T)J, 

AAK) - J g w _ ^ggp^ ^ T) - Bf (A e , Me , D) ■ (6) 

Above, uf and rip are the Fcrmi-Dirac distribution functions for the electrons and positrons: 

n F (E e ,fi e ,T)= Ec _l c - and fi F (E e) fj, e ,T) = EeW (7) 

e t + 1 er+1 



and A e = y/p 2 + m 2 is the electron energy. Note that the denominators in Eq. © give rise to higher than first 
corrections in the fine structure constant a = e 2 /Air, as can be seen by the expansion 

1 1 , M* 2 ) 5 



(A-P) 2 -(A 2 ) 2 /4 (A-P) 2 \iK-P) 2 , 

= Jk A pJ 2+ ° {a2) + --- (8) 

The 0(a 2 ) term would allow the unphysical decay of the plasmon into e + e~ pair through the imaginary part of the 
pole. The prevention of undesired 7* — > e + e~ decays entails a calculation that includes higher order diagrams which 
serve to increase the e + e~ pair mass and block such unphysical decays as explained in Ref. |l4j| . Thus, to 0{a), it is 
necessary to drop all but the first term on the right hand side of Eq. (JHJ) . 
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B. Polarization Functions 



In order to study the individual contributions of the different polarizations, it is convenient to express the vector 
part of r M „ in terms of its longitudinal and transverse components. (Following Ref. ^4|, we wu l utilize the term 
"photon" for the transverse com pon ent and "plasmon" for the longitudinal component.) The convention for these 
operators is adopted from Refs. |20ll21|. The axial part must be treated separately as it does not lend itself to such 
a decomposition. Since the polarization tensor is proportional to V^„, i.e., 

V^(K) = ( ^Xll flv (K), (9) 

the scalar polarization functions G{ujt, k) and F(u>l, k) can be defined through the relation 

W V {K) = G(cj T ,k)P» v + F(cj L ,k)P£ v , (10) 

where the tensor structure of IP 1 ' is captured in the polarization tensors P£ u and P^ v . The components of the 
transverse and longitudinal polarization tensors are 

„„ f for \i or v = 

Pt = i,j = 1,2,3 



L — 



K 2 



The polarization tensors satisfy the properties 

Pt^ p Pl P v = , P^ p Ptpv = ~Ptv 

P£ p Pl p „ = -PL , P£ M = -l , 

For completeness, we note that P pv can also be written in the form of a tensor product |l4j : 

1 
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Utilizing these properties, we have the relations 
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F(co L ,k) = P^U liU {K) = ^-U 00 (K). 
Performing the angular integrations in Eqs. JSJ) and <|15l) . the scalar functions G and F take the form 
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(15) 



(16) 



where v = p/E e is the velocity of the electron. These polarization functions are compactly defined as 

IL T (K) = G{uT,h), 

Xl L (K) = lf (K) = (k 2 /K 2 )F(Lu L ,k) (17) 

in Refs. [U[2l|. 

In order to calculate the axial part, it is convenient to write A^ v in terms of a simple tensor and a scalar function, 
as was done for Vp V previously. We would like to keep manifest Lorentz covariance in our expression, since frame 
fixing would require special care later when the emissivities are calculated in the frame of the heat bath. We perform 
calculations in the Lorentz gauge (9 • A = 0) which is suggested as a convenient choice by the equations of motion in 
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an effective theory of a massive spin-l particle (other choices j22| lead to the same result). From the tensor structure 
of in Eq. JfJJ), it is clear that only the transverse and longitudinal components of the integral 



d 3 p 

2E e (2ir) 3 (K-P) 2 



(n F - n F ) (18) 



contribute. The term "transverse" here refers to a 3-vector transverse vector, i.e., all three polarizations are 4-vector 
transverse to . As a result, only the longitudinal component survives the d 3 p integration (this is easy to see if we 
temporarily rotate the photon momentum along the z-axis). The integral above is then oc e' 3 ''', which in a general 
frame for a massive vector particle of mass M has the form j2^| 

XM'kMJ' { ' 

Using this property, it is clear that only the transverse polarizations contribute to the cmissivity. Thus A^ v can be 
recast as 

A »u = ^A^P^^g^^K^ (20) 



V2e Q 

where VK 2 = vTIt is the mass of the vector particle and the scalar function IT^ (K) was chosen to coincide with that 
given in Ref Hj: 

rr 2aR2 I j P 2 \^t, ( Lur + vk\ uj^-k 2 i, _ . 

UAiK) = J d ^[^ 1 H^^J"^^J (nF_nF) - (2 ) 

C. Squared Matrix Elements 

The squared matrix element for neutrino pair production in a plasma can be computed using 

3 

\ M \ 2 = J2 e^e^r'^r^L^ , (22) 

A=l 

where the summation is over photon polarizations A and the outgoing lepton tensor for neutrinos is 

L a = 8[<3i Q Q 2/ 3 + Q\pQlu — 9ap(Ql ■ Q2) + itaPpaQiQZ] ■ (23) 

The sum over polarizations can be performed using 

A=l 

From the symmetry properties of T^ v = + iA^ v and L^ v — + ia^ v ', it is easy to deduce the combinations that 
contribute to \M\ 2 . Explicitly, 

\M\ 2 = -V^V^v a p-A^A^v a p^{y^A^ a -V^A^)a a p, (24) 

where the transversality of r Miy , namely that K^T ^ v = 0, was used. The first term can be further decomposed into 
its transverse and longitudinal components to yield \M.\ T and l-M^. The second term is purely axial and will be 
denoted by while the last term is the mixed vector-axial term It will be shown later that this last term 

is purely transverse-axial (i.e. without longitudinal contributions) and does not contribute to the total emissivity and 
rate, but only to the differential cmissivity and rate. 

The squared matrix elements from the transverse, longitudinal, axial, and mixed vector-axial channels are: 
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(25) 



where the symbol (■ • • ) indicates that the appropriate spin averages and summations have been performed. These 
expressions are central to the differential and total emissivities presented in this paper. 

The main contribution to the total emissivity comes from the transverse mode, the longitudinal component becoming 
comparable to the transverse component only in the extremely degenerate limit fbH Il6l | . Contributions from the axial 
channel are orders of magnitude smaller than the transverse and longitudinal parts. The mixed vector-axial total 
emissivity vanishes completely. For all practical purposes, therefore, it is sufficient to consider only the transverse 
and longitudinal parts, while the axial and mixed vector-axial parts can be safely neglected. 



III. DIFFERENTIAL AND TOTAL EMISSIVITIES 

The total emissivity or the total energy carried away by the neutrino pair per unit volume per unit time can be 
computed from 
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where hb(uj,T) is the Bose-Einstein distribution function for photons or plasmons and the sum is over appropriate 
polarizations. The factor Z(k) arises from the residue of the pole in the propagator: 
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In vacuum, the free field solution has f(k) = m and therefore the denominator in the above expression is 2fco(k) = 
2w(k). In medium, various conventions for this correction can be found in the literature [III23. Here, Z(k) is chosen 
such that 
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Setting f(u>, k) = Ht(lo, k) in the case of transverse polarization results in 



Z T (k) 



<9Hr(wT, k) 



dwj, 



(29) 



A straightforward calculation then gives 
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where v = p/E e . Our choice 
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2u>Lvk \&l ~ vk j 
is slightly different from that used by Braaten 01 ■ Specifically, 

Z L = J*/(u>l - k 2 ) x Zl s . (32) 
This difference results in some minor differences in the expressions associated with the longitudinal emissivity. 



A. Transverse Differential Emissivity 

Integrating Eq. i|26[) over the 3-momentum of the photon, the transverse differential emissivity is given by 
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where 9 is the angle between the two outgoing neutrinos of energy E\ and E 2 . The dispersion relations are 

U) T {k) = k 2 +n T (io T {k),k) 

k = \Je 2 + E 2 +2E 1 E 2 cos9. (34) 

Equation Ij33() is the basis upon which suitable forms of the differential transverse emissivities can be obtained, 
since the energy <5-function can be used to advantage. For example, integrating over one of the outgoing neutrino 
energies (say, E 2 ) yields the differential emissivity as a function of the other neutrino energy [E\ = E) and the angle 
9 between the pair neutrinos. Explicitly, 
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where Jt(E,9) is the Jacobian of transformation from dE 2 to d[E 2 — lot(E\, E 2 )) originating from the (5-function 
integration: 
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For computational ease, it is advantageous to cast the dispersion relation in the form 
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which can be used to find lot(E,9) iteratively with Ht ~ lu 2 (see Eq. (|59|l 'l used as a first guess. In addition, the 
transcendental dispersion relation in Eq. (|34(l has to be solved iteratively for every choice of E and cos#. Using 
Eqs. lfT?)|> and IpTfjl . an explicit expression for Jt(E 7 9) can be derived. The result is 
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Alternatively, the energy <5-function in Eq. (|33|l can also be used to eliminate the angle between the two outgoing 
neutrinos. This leads to 
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where Jt(Ei, E 2 ) is the Jacobian resulting from the integration procedure. Its explicit form is given by 
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B. Longitudinal Differential Emissivity 



The procedure outlined for the calculation of the transverse differential emissivity can also be used to compute the 
longitudinal differential emissivity. However, there exists a difference between the two cases. For the longitudinal 
component, k < k max , so that the condition E\ + E 2 < k rnax has to be satisfied. This can be ensured by using the 
multiplicative factor 0(/C = k max — E\ — E 2 ). Utilizing the energy S- function to perform the integration over the 
Euler angles, we obtain 
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We first integrate over one of the neutrino energies to obtain the differential emissivity depending on E and 9: 
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On the other hand, performing the angular integration analogously to the transverse case gives 
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where cos6» = (k 2 — E\ - E 2 )/(2E 1 E 2 ). The factor qUe^ - (w| - fc 2 )) in Eq. gSJ) accounts for the fact that 
| cos 6*| is bounded by unity. 

C. Axial Differential Emissivity 

Inserting (|A^a| 2 ) from Eq. (|25|) into Eq. 1)26(1 and integrating over the photon momentum, the axial contribution 
to the differential emissivity is given by 
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Integrating over one of the neutrino energies yields 
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(2E - lo t ) 2 + E{uj t - E)(l + cos 



E 2 + (lu t - E) 2 + 2E{u T - E) cos ( 



where the residue factor Zt and the Jacobian Jt(E, 9) are given by Eqs. (|3(J[I and (|38|) . 

The corresponding expression for the differential emissivity as a function of the two outgoing neutrino energies is 



2 



d 2 Q A G ^f {C ^\z T{k)nB{ ^ T )^ T ,k)E 1 E 2 



dEidE 2 16n 4 a 



(k 2 +E 2 -E 2 )(k 2 -E 2 +E 2 ) 

ElE2 W 2 



x J T (E U E 2 ) e(4S 1J B 2 - IT T ) , (49) 

where Jt(E\, E 2 ) and Q(AE\E 2 — Ilx) originate from integrating over the <5-function, and Jt{E\, E 2 ) is given in Eq. 
60J|. 

D. Differential Emissivity from the Mixed Vector-Axial Channel 

The differential emissivities in this case are 

d3 QM G 2 F J2 f C f A C i (]A , n A (o; T ,fc)n r (^ r ,fc) 22f , 
d El dE 2 dcos9 ~ 8^ Z T {k)n B {u T ,T) E X E 2 {E 2 -E 1 ) 

x (1 - cos 6) S(uj t -Ex- E 2 ) , (50) 

d 2 Q M G 2 F J2fG A C v U A (ui T ,k)n T (uJT,k) 2 2 

— = Z T (k)n B (uj T ,T) E (w T -E) {u T -2E) 

dEd cos 9 87r 4 a k 

x (l-cos0) J T (E,9), (51) 
<?Qm_ _ Gl E f Cjcjr U A (. T , k )IlUuT, k) ^ _ 



dEidE 2 167r 4 a 

x MEi, E 2 )e{iE 1 E 2 -U T ). (52) 

Integrating over (E, 9) or {E\, E 2 ) confirms that the total emissivity Qm — 0, a result that was anticipated from the 
V\ <-» v 2 antisymmetry of the squared matrix element in Eq. 12oli . 

The differential emissivities for the transverse component in Eqs. I|35fl an d 139|) . the longitudinal component in 
Eqs. I|43|) and l|45|). the axial part in Eqs. H48fl an d l|49f) . and the mixed axial- vector part in Eqs. 151|) and l|52|) are 
among the principal new results of this work. 
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E. Total Emissivity 

The total emissivity can be obtained from Eq. I|26|) by performing integration over the entire phase space and 
summing over all three outgoing v flavors. Accounting for a factor of 2 for the two transverse photon polarizations, 
the transverse emissivity takes the form 



x IL t (u>t, k) 



d3k z T (k) Al 



d 3 



(12 



2w t (2tt) 3 v ' 2£ 1 (2tt) 3 2E 2 {2tt) 3 
(k-qi)(k-q 2 ) 



{Ei + E 2 )n B (uT, T)(2ir) 4 5 4 (K - Q 1 - Q 2 ) 



E1E2 



k' 2 



2 v'i 



dk k Z T (k) n B (u T , T) (u T - k ) 
Similarly, we obtain the longitudinal emissivity for single plasmon polarization as 



!L = 



967r 4 a 



dk k 2 Z L {k)n B (to L ,T) (ul-k 



2>:s 



(53) 



(54) 



where k max is the light-cone limit (i.e., the maximum k for which k < ujl) for the longitudinal plasmon. The results 
in Eqs. I|53|) and l|54|l agree with those of Ref. 01 modulo slight differences in the conventions used for Z^. 

Integration over the two outgoing neutrino momenta and angles in Eqs. (|25|l and l|26|) yields the total axial emissivity 



-2 



Q, 



r>i \r- /-if 

U F 2^f °A 

487r 4 a 



dkk 2 Z T (k) n B {uj T ,T) U 2 A (uj T ,k) (w| - k 2 



(55) 



which agrees with the result given previously in Ref. jl4j . As noted earlier, the total mixed vector-axial emissivity 
Qm = 0. 

A cross-check of Qt, Ql and Qa is afforded by integrations of their respective differential emissivities and by the 
results reported in Ref. pL4j | . We present such checks in the following section. 



IV. RESULTS AND DISCUSSION 



For the most part, we will present our results for the total and differential emissivities as a function of the mass 
density of protons in the plasma, p B Y e = m p n e , where m p is the proton mass, Y e = n e /n B is the net electron fraction 
(n B is the baryon number density) , and n e is net electron number density 

1 f°° 

n e (T, /i e ) = — dpp 2 (n F -n F ), (56) 
Jq 

which is simply the difference between the e~ and e + number densities. (The quantity p B Y e is the same as pj fi e used 
in prior works including Ref. [13.) Given n e , this expression can be inverted to determine the chemical potential 
p e at a given temperature. In connection with supernova simulations, however, it will be more natural to examine 
emissivities as a function of the baryon mass density p B = m p n B . 

The differential and total emissivities require as inputs the transverse and longitudinal dispersion relations u>t,l {k) 
in Eqs. (|34() and (|42() , which are transcendental relations that involve the polarization functions in Eq. I|17fl which 
are in turn expressed as integrals that involve uix,L(k). In addition, the residue factors Zt,l hi Eqs- l|30|) and 1311) 
also require tUT,L{k) as inputs. 

With the help of the approximate expressions developed earlier in the literature (see, for example, [hH Il6l | and 
Sec. II V A1 below) . we have calculated the exact dispersion relations numerically by using iterative techniques. En- 
ergy integrations were performed by using standard Newton-Cotes algorithms. Here an appropriate multiple of the 
temperature was used as the high energy cutoff. Angular integrations were performed by employing Gauss-Lcgendrc 
quadrature with N = 16, 32 and 64 points, all of which converged to the same result. 

Figure [2] shows the exact results of transverse and longitudinal dispersion relations for temperatures T = 10 8 K 
and T = 10 11 K, respectively. These calculations are numerically cumbersome and time consuming, particularly in 
the case that differential emissivities are needed. We therefore turn to assess the extent to which the approximations 
developed in Ref. fl4T | reproduce the exact results over a wide range of temperature and density. 
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FIG. 2: Top panel: Dispersion relations in Eqs. (1341 calculated numerically using the exact polarization functions in Eqs. 1171 
and 1421 1 at the indicated temperatures and densities. The middle and bottom panels show the relative accuracies of the 
longitudinal and transverse dispersion relations calculated using the approximate polarization functions in Sec. IIV Al 



A. Approximations for Z, II, J and j3 

For completeness, we collect here the approximate relations developed in Ref. 01 . These relations are particularly 
helpful in developing new expressions that are needed in the calculation of differential emissivities. For the transverse 
mode, the residue factor Zt and polarization function T\t are: 



Z T (k) 



2lo 2 t {lo 2 t -v 2 k 2 ) 



3L0 2 LO 2 + ( w 3 + fe2)( w 2 _ „2 fc 2) _ 2W 2 ( w 2 _ fc 2) 



(57) 



where 



n T (w T , k) = uItt? 



vlk 2 LO T 



In 



( wr -f 




\lO T - 


v*k J 


+ n F ) , 





dp k(}- v - ]{nF 

4a [°° p 2 (hv 2 
-J a dP K{—~ 

LJl/u p . 



(58) 



(59) 

(60) 
(61) 



In Eq. 1)59(1 . lo p is the plasma frequency - the lower limit for the value of plasmon mass (in this case transverse, since 
IIt > Wp). The quantities wi and are auxiliary variables that render the expressions for Zy and Ily (V = T or L) 
compact. 

Equation (|58^l can also be used to obtain approximate results for Eqs. (|3*%jl and l|4()|l. For Jt{E, 9), we get 



Jr(E,t 



lo t (1-^ 9) ) 



with 



lo t {^ - 1)4 B ' S) + E(l - cos0)(l - ^# e) ) 

( LOT + V* k 



3 _ 



2v^kLOT 



111 



(62) 



(63) 
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and for Jt(E±, E 2 ), we obtain 



Jt(E 1 ,E 2 ) = 



E\Ei 



1 



1 - 



k 2 



with 



9 < 



4 vlk 2 



1 / u* k 3wt \ / + v * fc 



6 V wt v*k I \ljt — v*k 



(64) 



(65) 



where = E± + E 2 . 

For the longitudinal mode, approximate expressions for Zl and II £ are: 



Z L (k) 
U L (u L ,k) 



LO 2 - k 2 3^2 - ( w a _ ^2) 



2v*/c 



hi 



oj l - v*k 



- 1 



(66) 
(67) 



Compared to Ref. ^4[, the presence of the additional factor Lo 2 L j{uj 2 L — k 2 ) in Zj^ stems from the differences in the 
convention adopted for Z^. In this case, the Jacobians are: 



and 



Jl{E,9) = 



Jl(Ei,E 2 ) 



p(E u E 2 ) 



kHi-pr ] ) 



k 2 (l- + lo l {u l -E(l- cose))^ 9) 



2vl 



2fc3" 



111 



ujl + v*k 
ijJl — v*k 



p[E u E 2 ) 



k 2 u 2 L — v 2 k 2 



2v* 
I 2 " 



E 1 E 2 uj l (3 {EuE2 



3<4 
2vl 



2fc3" 



111 



OJL + V*k 



v 2 k 2 



2v* 



with ujl = Ei + E 2 . Ref. [l4T | provides an approximate expression for Ua(u>t, k): 

' - k 2 2>ul - 2(cj| - k 2 ) 



U A (u> T ,k) = ui A k- 



•j 2 k 2 



with 



2a [°° p 2 2 2 _ 
-J *^(1-3«)(»f-»f)- 



(68) 
(69) 

(70) 
(71) 

(72) 
(73) 



We turn now to numerical results employing these approximations. The lower panels in Fig. [21 show the relative 
difference between results calculated using these approximations and the exact, but numerical, calculations of the 
dispersion relations. At the temperatures considered, the difference between the two dispersion curves is at most 
of order of 10 -5 for the transverse case and decreases further as k increases. For the longitudinal case, the relative 
difference is somewhat larger, but still negligible. These approximations greatly accelerate computations without 
significant loss of accuracy. 



B. Checks of Differential Emissivities 



The expressions for the differential emissivities in Sec. IIIII can be integrated over (E 7 8) or (Ei,E 2 ), in order to 
recover results for the total emissivities calculated independently from Eqs. I|53|) and l|54|) . respectively. This is not 
only a useful check of results derived for the differential emissivities, but also allows for a qualitative understanding of 
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FIG. 3: Total transverse emissivities Qt from Eqs. 153H versus baryon density at the indicated temperatures. The symbols 

'x' and '+' show results of J d ^ d ^ s e dEd cos 9 and f d d E ^] S2 dE\dEi using the differential emissivities from Eqs. 1351 1 and 1391 ■ 
respectively. The vertical dotted lines mark the turn-on densities at which Qt abruptly changes its density dependence. 

photon and plasmon decay under various physical conditions of the plasma. In this section, we analyze the transverse 
cmissivity in some detail. The longitudinal and axial emissivities are discussed in Sec. IIV C 31 

Integration over energies and angles were performed by substituting the approximate expressions for IFr, Zt, and 
J T (E,9) and J T {E 1 ,E 2 ) in Eqs. (JSTJl, ||5HJ, (EH, and $5^. The results are shown in Fig. |3J It is gratifying that 
there is excellent agreement between the emissivities calculated in two different ways. In addition, they also agree 
with those published in the literature 0, 0, . 

While the transverse total emissivities shown in Fig. |3| increase rapidly with temperature, their behavior as a 
function of density is more complex. At a fixed temperature, the basic characteristics to note are: 

(1) Qt is independent of the density psY e until a turn-on density (pBY e )to is reached, 

(2) For densities larger than this turn-on density, Qt exhibits a power-law rise until a maximum is reached at 
(pB^ e )pcak, and 

(3) For pbY c ^ {pBYe)pcak, the fall-off with density is exponential. 

In the next section, we proffer both qualitative and quantitative analyses of these features. 

C. Qualitative Behaviors 

In order to gain a qualitative understanding of the basic features of Qt in terms of the intrinsic properties of the 
plasma, it is instructive to inspect the behaviors of the chemical potential fi e (together with T, this determines n e ) 
and plasma frequency lv p (this is the characteristic energy scale generated by interactions in the medium) as T and 
psY e are varied. Fig. shows p e (lower panel) obtained by inverting Eq. I|5t)[) and w p (upper panel) calculated from 
Eq. (|59|l . The results in Fig. [31 are readily interpreted on the basis of the trends observed in Fig. 21 To establish 
the main points, we focus on the transverse emissivity Qt in this section. The sub-leading longitudinal and axial 
contributions can be understood in a similar fashion (see Sec. IIV C 3|) . 

Noteworthy features of the plasma frequency in panel (a) of Fig. 01 are: 

(1) ujp is independent of psY e till the turn-on density (pBY e )to is reached (this is at the root of why Qt is constant 
for p B Y e < (p B Ye)to), and 

(2) u)p shows a power-law increase for psY e > (psi^)to, the index depending both on the extent to which the 
plasma is in the non-degenerate, partially degenerate or degenerate regime and on whether electrons are relativistic 
or nonrelativistic. 

The bold and light portions of the various curves in panel (b) of Fig. 0] mark the regions of densities for which 
p, e > T and \i e < T) , respectively. Inasmuch as p, e ~ T indicates partial degeneracy of the plasma, the bold and light 
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FIG. 4: The plasma frequency co p (upper panel) and the chemical potential /i e (lower panel) as functions of density for select 
temperatures. The vertical dashed lines in panel (a) indicate the densities at which u p = T. The vertical dotted lines mark 
the turn-on densities at which the plasma frequency and the neutrino emissivity (Fig. [3J abruptly change their behavior from 
being independent of density. 

portions refer to the degenerate and non-degenerate conditions, respectively. For reference, the electron mass, which 
when compared with /i e or T determines the degree of relativity, is marked by the horizontal dotted line in this figure. 
The vertical dotted lines show the respective locations of the turn-on densities. 

In order to assess the role of positrons in determining the behavior of the plasma frequency, we show the electron 
and positron fractions versus density in Fig. [BJ For a given temperature, the positron fraction begins to decrease with 
increasing density or chemical potential and vanishes exponentially in the limit of complete degeneracy. It is interesting 
that the positron fraction is of order 15% at the turn-on densities for each temperature shown in Fig. [3J Although 
this is small compared to the electron fractions at these densities, it will be necessary to account for the presence 
of positrons to gain a quantitative understanding of uj p as a function of density and temperature. Since varying 
degrees of degeneracy and relativity are encountered in the wide ranges of density and temperature considered here, 
we consider below those limiting cases that help us to understand the main features of oj p and Qt- 
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FIG. 5: Electron and positron fractions in the plasma versus density for select temperatures. The vertical dotted lines mark the 
densities at which the plasma frequency (Fig.^J and the neutrino emissivity (Fig.[3J abruptly change their density dependence. 



1. Plasma Frequency 

We begin by enquiring why ui p is constant with increasing psY e at low densities. Consider first the cases of 
T = 10 11 K and 10 10 K, both of which exceed the electron mass (m e = 5.92 x 10 9 K). The plasma is thus in the 
relativistic regime in which the net electron density and the plasma frequency are well approximated by 

A*e ( 2 , rrn2 2 \ PeT 2 
1 IT 



3tt 2 V / 3 ' 

2 4a/ 2 7r 2 T 2 \ 4naT 2 

< = + — - (74) 

In writing the rightmost relations above, we have used the fact that p e <C T {i.e., the plasma is non-degenerate) as 
is evident from the lower two curves in panel (b) of Fig. The plasma frequency is therefore effectively independent 
of psYe and is given by 

(w p )pi~r/10, (75) 

where the subscript "pi" denotes the plateau value. With increasing density, the chemical potential p, e rises linearly 
with density (see Eq. J?1J) until the plasma enters the partially degenerate regime in which /z e becomes comparable 
to T. In this regime, lo p abruptly changes its behavior from a constant to a power-law behavior. An estimate of the 
density, (pBY e )to, at which this turn on occurs can be obtained by setting p e ~ T in Eq. I|74(l . whence we have 

(n e ) to ^ tt^ or 1 

3 } Relativistic (76) 



( PB Y e ) to ~ 8 x 10 7 (j^j g cm- 



For T = 10 n K (8.62 MeV) and 10 10 K (0.862 MeV), Eq. (J75J) yields log 10 (p B Y e ) to = 10.7 and 7.7, respectively, in 
agreement with the results shown in Fig. 01 

Further increase in density results in the plasma becoming degenerate as p e 3> T. Using the Sommerfcld expansion, 
we find the net electron density and plasma frequency in this region to be 



= U + ^-(2p% + rnl) (77) 



2 4a/ 2 k 2 T 2 



3m 2 „ 2 
e - + 2p F v F 



^P>f . (78) 
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where Ep — (p 2 F + m 2 ) 1 / 2 and vp = pp/Ep. The rightmost relations here refer to the case in which the thermal 
contributions are small, i.e., T -C p e — Ep. In this regime, the plasma frequency varies with the net electron density 

1/3 . 1 

as uj p oc n e , since vf — > 1. 

For T = 10 9 K, the region of densities in which u p is constant lies in the classical regi me, in which T <C m e and 
me — Me ^ T. Including the contributions from positrons (this was neglected in Ref. [22|), we find that the net 
electron density and the plasma frequency in this non-degenerate and non-relativistic limit are 



n, 



U) . = ^ mh (ti)('JliL) : " 2 e -'./T( l _ 5T 



Anan e ( u P \ ( 5T 
coth 



For low densities at T = 10 9 K, p e <C T which renders u) p to be independent of p e and to be a function of T alone and 
thus independent of n e . At this temperature, Eq. (|80|l yields log 10 uj p ~ —2.95, in close agreement with the exact 
result shown in Fig. The turn-on density at which uj p begins its characteristic power-law rise with density can be 
obtained by setting fi e ~ T in Eq. I|79|) . At this temperature, we find log 10 {pBY e )t ~ 3.2, in excellent agreement 
with the exact results shown in Fig. 0] 

For pbY c > (pBY e )to, we observe that uj p oc (psY^) 1 ^ 2 as the plasma begins to enter the partially degenerate regime 
from the non-degenerate regime. For densities characterized by /i e 3> T, the plasma frequency exhibits the same 
degenerate behaviour as for the higher temperature cases, i.e., lj p oc (p^Ye) 1 / 3 . 

For T = 10 8 K, the region of densities in which oj p is constant with density lies well below the lowest density shown in 
Fig. 0] Note that the (ps^e) 1 / 2 rise, typical of a non-degenerate plasma, is also observed for a degenerate (p e >> T), 
but semi-relativistic (m e ~ T) plasma. The feature that p e is nearly constant until a density of 10 6 g cm -3 is 
connected with the fact that in this highly degenerate case, 

m 2 + (3^ 2 n e ) 2 / 3 ] 1/2 (81) 

begins to rise significantly with density only for densities that exceed the threshold density of m 3 /(37r 2 ). This situation 
does not occur for the other temperatures shown in Fig. 01 to the degree that it does for T = 10 8 K. 

We turn now to determine the density at which tu p = T . Consider first the degenerate case in which fi e 3> T and 
the plasma is rclativistic (vp — > 1 and /i e — * Pf)- Ignoring the small temperature dependent corrections, Eqs. I|77f) 
and l(7%|) can be combined to yield 



Rclativistic (82) 



Results from this expression agree closely with the numerically calculated exact results (see the vertical dashed lines 
in panel (a) of Fig. |2J for all but the lowest temperature of T = 10 8 K for which Eqs. 1)77(1 and 1(78(1 can still be used 
to advantage, but with vf = pf/Ef- The leading order solution of the resulting cubic equation in p 2 F leads to 



I v \ m p m e T 2 
(PBY e ) u _ T = 



Aira 



~ 1.22 x 10 9 (-^—] g cm" 3 . 
which agrees with the exact result at T = 10 8 K. 



Nonrelativistic (83) 



2. Transverse Emissivity 

The results of the previous section enable us to organize the discussion of the basic qualitative features of the 
transverse emissivity Qt, which is largely governed by the temperature T and the plasma frequency uj p (fi e ,T). 
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Depending on the extent to which the plasma is degenerate and rclativistic, uj p is either a simple function of p e or T 
alone, or an involved combination of both p e and T. The density regions where u p < T and uj p > T identified in the 
previous section are particularly useful here as they enable us to examine the limiting cases u p <cT and u> p T in 
which Qt in Eq. (|53|> can be cast in physically transparent forms [l4j . 

In the case of T > w p , the main contribution to Qt comes from high photon momenta. To a very good approxi- 
mation [T^j . 

where £(3) ~ 1.202 is Ricmann's Zcta function and tot is the transverse photon mass (this can be read off from the 
high energy limit of the transverse polarization function, i,e., ILt(u>t, k) k ^°° > m T ; the transverse mass lies in the 
range u p < tot < \/3/2 Up). With tot = u p , 

(Q T )pi ~ 2.8 x 10 24 u 6 p T 3 erg s~ : cm" 3 , (85) 

where lo p and T are in units of MeV. As discussed in the above section, ui p is independent of density for psY e < 
(pBYe)to- Consequently, Qt is also density independent up to (pBY e )to, a trend which is maintained at all temperatures 
whenever u p is density independent. Use of Eq. (|85|l with uj p from Eqs. (|74|l and 1)80 ll yields excellent agreement with 
the plateau results at all temperatures shown in Fig. [3] 
For T <C u p , Qt takes the form Q 

This expression allows us to qualitatively understand the power-law rise of Qt for densities pBYe > (ps^e)to- In the 
classical or non-degenerate regime, lo p cx nj 2 , whereas in the degenerate regime, ui p cx nl^ 3 . As a result, 

Qt oc (p B r e f , 2.5 < p < 3.75. (87) 

For all temperatures considered, the peak value of Qt occurs in the relativistic and degenerate regime in which 
T <C Up. Noting that at the peak w p = 7.5 T, we can utilize Eq. (|51)|> to obtain 

/ t \ 9 



(QT)pcak ^ 2.9 x 10^ (^^J crgs^ 1 cm" 3 , (88) 

which yields the peak values in Fig. |3| The density at which the peak occurs may be obtained by inverting Eq. I|78|) 
with the result 

1 /3r^ 3/2 



Using the fact that lu p = 7.5 T at the peak, we arrive at 



(p B Y e ) 



3375 / 3 \ 1/2 / T \ 3 
pcak _ 32^ J mp VMcV J 

- 1-8 x 10 13 (^) gem" 3 , (90) 



which accounts for the results at all temperatures shown in Fig. 

As the net electron density increases to values well in excess of (pBY e ) pea k, the plasma frequency begins to become 
significantly larger than the temperature. In this case, the emissivity is exponentially damped by the Boltzmann 
factor e~ u}p l T as seen in Fig. |3 



3. Comparison of Transverse, Longitudinal, and Axial Emissivities 



The qualitative analysis of the longitudinal and axial emissivities can be carried out along the same lines as that 
for the transverse emissivity. Figure shows the individual contributions from the three channels at T — 10 n K and 
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T = 10 9 K, respectively. The symbols "x" and "+" show results obtained by integrations of the differential emissivities 
in Sec. ED 

The plateau value of Ql can be inferred by considering the limit T 3> lo p . In this case, 

QL ^f/% v.Kr, (9i) 

which shows that Ql is suppressed relative to Qt by a factor uj 2 /T 2 [l^. In the relativistic limit, the coefficient 
A(v* -> 1) ~ 0.349 Using the result uj 2 p = AnaT 2 /9 from Eq. JJJ) that is valid for /i e < T, we find 

-.15 // ^ A „_„ „„,-3 „-l 



(Q L ) pl ~2.14x 10 15 J crgcm^s" 1 , (92) 

which is in fair agreement with the exact result for T = 10 n K shown in Fig. 

The turn-on density and power-law rise of Ql with density are similar those of Qt- As for Qt, the peak value of 
Ql occurs in the regime T « w p . In this case, Ql has nearly the same functional form as Qt in Eq. I|86|) . but with a 
slightly different numerical coefficient (for vp — > 1, Ql — V% Qt) EH- For a fixed temperature, therefore, the peaks 
of Ql and Q T occur at the same density with nearly the same values. 

We discuss now the relatively small contribution of the axial channel and focus on its peak value, since the plateau 
region for this case is at densities substantially lower than those of interest here. As shown in Ref . 01 > 

Qa * nci 3V2^ ( 1 + U ) ^iT^e-^ (93) 



96?r 4 a \ 5 

for T < w p . In the degenerate and relativistic limit, (T <C /i e and —* vp — > 1), it is straightforward to establish 
that the factor uja — (a/S^) 1 / 2 ^. This implies that the peak of Qa occurs at oj p = 6.5 T. Exploiting this result, we 
arrive at 

23 f T \ „,.„ „,„-3 „-i 



(QA)pcak ^ 5 x 10^ 1 — 1 erg cm^ s" 1 , (94) 

which nicely reproduces the exact results at T = 10 n K and T = 10 9 K shown in Fig. EI 

Eqs. H77f) and J7SJ, coupled with ui p = 6.5 T, yield the density at which Qa peaks. Explicitly, 

(PBY e ) pcak ^lA7xlO^(J-) gem- 3 . (95) 



This density is slightly smaller than the density at which Ql and Qt attain their peaks, since u p = 6.5 T at the peak 
of Qa, while tu p = 7.5 T for the other two channels. 

Our results for Qa differ by orders of magnitude with those in Fig. 2 of Ref. 0| , although the formal expressions in 
both works are identical. For example, at the maximum density of log 10 ps^e = 14 at T = 10 11 K shown in Ref. [l4j . 
log 10 Qa — 32.8 whereas we obtain ~ 31. At the lowest density of log 10 pbYe — 10.5 displayed there, log ln Qa — 24, 
while we find ~ 23.4. This latter result can also be verified from the analytic expression in Eq. (34) of Ref. [lj], which 
is valid in the regime T 3> w p . We have been unable to resolve these large discrcpcncics in the numerical results for 
Qa- Since we agree with all of the formal results derived in Ref. 0|, we attribute these differences to a possible 
numerical error in the calculations of Qa in Ref. |14| . We wish to emphasize, however, that these discrepencies are not 
of much significance, since Qa is orders of magnitude lower than both Ql and Qt for the densities and temperatures 
of interest here. 



D. Features of Differential Emissivities 



An interesting dependence of the differential emissivity on the angle between the two outgoing particles is worth 
noting. If, for example, the differential emissivity from Eq. (|35|1 is integrated over the outgoing energy, the only 
remaining dependence is on the outgoing angle 9. For a massive particle decaying into two particles, we would expect 
"back-to-back" decay, i.e., the two outgoing particles prefer to be emitted at 8 = 180°, for which cLQt /d(cos 9) should 
exhibit a maximum. The massiveness of the photon can be gauged by the photon effective mass which strongly 
depends on the properties of the medium, such as its temperature, chemical potential, and net electron density. 
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FIG. 6: Individual contributions from the transverse, longitudinal, and axial channels to the neutrino emissivity. 




FIG. 7: Angular dependence of the differential emissivity for two extreme cases. Panel (a) shows the dominance of energy 
emitted at a small outgoing angle between the neutrinos for densities in the plateau range psY e = 10 — 10 10 g/cm 3 , while panel 
(b) depicts a situation in which "back to back" emission is preferred for densities in the peak range psY e = 10 12 — 10 13 g/cm 3 . 



Effectively, the photon mass is roughly equal to oj p and its kinetic energy is of order T. For uj p much larger than the 
kinetic energy, we can expect aspects associated with the decay of a massive particle. However, if its mass is negligible 
compared to the kinetic energy, we should see a significantly different behavior. 

These different behaviors arc illustrated in the results presented in Fig. [7| In panel (a), the "light" photon (i.e. 
ojp <C T) decay is shown. In this case, the main contribution to the total emissivity comes from neutrino pairs with 
small outgoing angles between them. This corresponds to a high Lorentz boost of the photon with respect to the rest 
frame where calculations were performed and the angle that was « 180° in the photon rest frame gets significantly 
shifted. When the photon effective mass becomes the dominant energy scale (u> p T), the result is the massive 
particle decay depicted in panel (b). 
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V. KERNELS FOR NEUTRINO TRANSPORT CALCULATIONS 



In this section, we discuss briefly how the differential emissivities calculated in this work enter in calculations of 
neutrino transport. The evolution of the neutrino distribution function / is generally described by the Boltzmann 
transport equation in conjunction with hydrodynamical equations of motion together with baryon and lepton number 
conservation equations. Ignoring general relativistic effects for simplicity (see, for example, Ref. |l2j| for full details), 
the Boltzmann equation is 

^ + v % + = B EA(f) + B NES (f) + b umU) + B T p(f) , (96) 

where F l is the force acting on the particle. In the source term on the right hand side, Bea(/) incorporates neutrino 
emission and absorption processes, Bnes(I) accounts for the neutrino-electron scattering process, B u j^(f) includes 
scattering of neutrinos off nucleons and nuclei, and Bxp(f) considers the thermal production and absorption of 
neutrino-antineutrino pairs. 

In discussing the contribution of the plasma process to Btp, we follow Ref. [l2T | in which neutrino pair production 
from e + e~ annihilation was considered in detail. Suppressing the dependencies on (r, t) for notational simplicity, the 
source term for the plasma process can be written as 



B(f{jii,E 1 )) 



l-/(A*i,£i] 



poo pi pZ-n / 

J El dE 2 J dfi 2 J d<t> 2 R p \Ej. , E 2 , cos 9 



l-fU 2 ,E 2 



(2tt) 3 

" /(Ml ' Sl) (2^/ EldE 2 J^ 2 ^ d^ 2 R a ^E 1 ,E 2 ,co S e^j f(fi 2 ,E 2 ), (97) 



where the first and the second terms correspond to the source (neutrino gain) and sink (neutrino loss) terms, respec- 
tively. Angular variables jjn = cos 6i and <pi arc defined with respect to the z-axis that is locally set parallel to the 
outgoing radial vector r. The angle 9 between the neutrino and antincutrino pair is related to 9\ and 6 2 through 



cos 6» = piiM2 + \J{1 - - cos(0!-0 2 )- 

The production and absorption kernels are given by 

R*(Ex,E 2 ,cos6) = / 7^ Z Y (k) T M 1 S 4 (K - Q 1 - Q 2 )(2^(\M | 2 ) , (99) 



(2tt) 3 w V 1 - n B (u,T)J %uE x E 2 

where the subscript Y stands for T- "transverse" or L- "longitudinal" . The factor £ accounts for the spin summation; 
£ = 2 for the transverse, axial and mixed cases, while for the longidudinal case £ = 1. 

The angular dependences in the kernels R" (Ex, E 2 , cos 9) are often expressed in terms of Legendre polynomials as 

R*(Ex,E 2 ,cos6) = f^^^^{Ex,E 2 )Pi(cose), (100) 
1=0 

p 

where the Legendre coefficients (Ex, E 2 ) depend exclusively on energies. 

From Eq. I|99(l . it is evident that the kernels are related to the neutrino rates and emissivities. We first consider 
the production kernel R p (Ex, E 2 , cos 9). The corresponding analysis for the absorption kernel R a (Ex, E 2 , cos 9) can 
be made along the same lines, but with the difference that ns is replaced by 1 — i%b- The neutrino production rate 
is given by 

r = c / wk ZAk) m m T > - - Q^(\Mf) 

d s qi d3< &. R P( EuE2}COse ^ (ioi) 



(2tt) 3 (2tt) 3 ~ 



which defines the kernel R p (Ex, E 2 , cos 9) and is to be identified with that in Eq. I|99|) . The emissivity Q can also be 
cast in terms of RP using 

Q = j ^-^(Ex + E 2 )RP (Ex, E 2 , cos 6). (102) 
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Equations (fTUT|) and (|H)2"|) can be inverted to obtain 

8tt 4 d 3 T 



R*>(Ex,E 2 ,cos6) - ^ ———— 

^ d3Q (103) 



E\E 2 (Ex + E 2 ) dE 1 dE 2 dcos9 ' 

Utilizing the results from Eqs. i|33|) . (|41|l . (|47|l . and (|50|l in Eq. (|1U3|> . the production kernels for the transverse, 
longitudinal, axial and mixed vector-axial parts arc: 



R p T (Ex,E 2 , cos 9) = 

R p L (Ex,E 2 , cos 6) = 

R P A (E 1 ,E 2 , cose) = 

R P M (E 1 ,E 2 , cos 6) = 



Z T (k) n B (ujT,T) 



G%J2f(Gy) 

2a lot 
x (l-cos0) 6(lv t - Ex - E 2 ) , 



IL T (m T , k ) {E 1 - E 2 ) 2 + E x E 2 ( 1 + cos ( 
E\ + El + 2E X E 2 cos 9 



(104) 



' Z L (k) n B ( UL ,T) K - fc 2 ) 2 ElE2 1 'g S % L -Ex- E 2 )Q(fC) , (105) 



2q 



Z T (fc) n B (uj T ,T) 



U 2 A (cu T ,k) (Ex-E 2 ) 2 + ExE 2 (l + cos( 



2a cjt 
x (l-cos^)5(a;r--E'i -£2), 



Z T (k)n B (oJT,T) 
x (1 - cos 6) S(uj t -Ex-E 2 ) 



LdTk 



E\ + E 2 + 2ExE 2 cos9 
(E 2 - Ex) 



(106) 



(107) 



A simplifying feature in these expressions is worth noting. Since the two outgoing neutrinos stem from the decay 
of a single photon or plasmon, an energy-conserving <5-function survives in the expressions for the kernels. This 

is fortunate, since it makes the calculation of the Legendre coefficients $j a (Ex,E 2 ) straightforward. The Legcndrc 
coefficients are determined from 



$f{Ex,E 2 ) 



d(cos 9)Pi (cos 9)R* (Ex , E 2 , cos t 



Integrating over the <5-function introduces a Jacobian and sets cos 9 to be 

cos9 



k 2 - E 2 - E 2 



2ExE 2 

As a result, the Legendre coefficients for the production kernels are: 



(108) 



(109) 



^(Ex,E 2 ) 



$ P l(Ei,E 2 ) 



^ P a(Ex,E 2 ) 



&m(Ei,E 2 ) 



ExE 2 



2a lvtExE 2 
x J T (Ex,E 2 )e(mE 2 ~U T )P l (cos9), 

2a M > By ^ ' uj l E 2 + E 2 +2ExE 2 cos( 

x J L (Ex,E 2 ) 0(/C) Q(4ExE 2 - (u>\ - k 2 )) P z (cos0) , 



(fc 2 + E 2 - E 2 )(k 2 - E{ + E 2 ) 



4k 2 



ZT(k)n ^ T , T) ^,k) 

2a lutExE 2 
x J T {Ex,E 2 )e{AE 1 E 2 -n T )P l (cos9) 7 



ExE 2 



(k 2 + E 2 - E 2 )(k 2 - E 2 + E 2 ) 



4fc 2 



/~i2 (~<f (~<S 



7 (1\ I ^ UA(0JT,k)U^(uJT,k) 

Z T [k) n B (uj T ,T) — ±-z (E 2 - Ex) 



2a ku>TExE 2 
x J T (Ex,E 2 )e(4ExE 2 -U T )Pi(cos6) 



(110) 



(111) 



(112) 



(113) 
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The corresponding absorption counterparts are 

E^(E 1 ,E 2 , cos 6) = 1 ~ n f"£ ) E? x (E 1 ,E 2 ,cose) (114) 

K(Ei,E 2 ) = 1 ~ n f^h l(E 1 ,E 2 ), (115) 

where x stands for T, L, A, or M. As for the differential emissivities in Sec. I1II1 the appropriate dispersion relations 
have to be solved iteratively for every choice of variables. 



Typical Neutrino Energies 



The ratio of the total emissivity to the total rate, Q/T, is a good measure of the mean neutrino plus anti-neutrino 
energy. Since the dominant contribution arises from the transverse photons, we focus on 



(E V u)t = -pr - - 



dk k 2 Z T (k) n B (oj T , T) (uj t - k 2 ) 



2\3 



k 2 

dk — Zr(k) ub(lot 7 T) (lo t 



(116) 



2\3 



k 2 ) 



results for which are shown by the solid curves in Fig. [S] The behavior of (E vV )t with density and temperature can 



be understood by examining limiting cases. Setting Zt 
for semi-quantitative estimates, we obtain 



1 and using the approximate relation lo t — k~ 



LO„ 



(E v „) 



T 



3 = 1 

OO 

£>iOW(iy) 



where y = lu p /T and K v is the modified Bcssel function of order v. 
estimates in the extreme relativistic and nonrclativistic cases: 



(117) 



Eq. (|117fl can be used to obtain numerical 



1.46 T 



for 
for 



(118) 



These results give a semi-quantitative account of the exact results shown in Fig. [S] The constancy of (E u „}t at low 
density and the power-law rise at high density arc chiefly due to a similar behavior of io v with density (see Fig. 

Since neutrinos stem from the decay of transverse photons in the plasma, additional insights into the magnitude 
of {E uS ,)t can be obtained by examining the mean transverse photon energies {lot) as functions of density and 
temperature. Explicitly, 



(lot) 



dk k 2 LOT ub{lot^T) 
dk k 2 nB(iOT, T) 



oo 

u p j=l 

8T 



K ±{jy) - K o{jy) 



3 = 1 



where the rightmost relation is obtained by using the approximate relation u> T — k 2 
relativistic and nonrclativistic cases, (lot) has the familiar forms 



(119) 



In the extreme 



(lo t ) 



2.7 T 



T 



for 
for 



T > L0 V 



T -C lo % 



(120) 



Results obtained by using the exact dispersion relation (the dashed curves in Fig. |SJ are matched well by these 
estimates. As expected, (£W)t ^ (wt) f° r T "C lo p . The suppression of (E vV ) T relative to (lot) for T » lo p at low 
density is principally due to the suppression inherent in the squared matrix element for neutrino pair production from 
the transverse photon in the plasma. 
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FIG. 8: Average transverse photon and neutrino energies (lot) and {E v „)t as functions of density and temperature. 



VI. COMPARISON WITH COMPETING PROCESSES 



In addition to the plasma process considered in this work, vv pairs are produced from the annihilation of e + e~ 
pairs (e + + e~ —►!/ + £/), the interaction of electrons with photons (e~ + 7 — > e~ + v + v), and the bremsstrahhmg 
process involving nucleons (n-\-n—tn + n-\-v + i>). In what follows, we wish to assess the relative importance of these 
processes in the ranges of density and temperature considered in this work. In view of the developments reported 
in Ref. \\M a nd this work, our results represent an update of similar comparisons made earlier in the literature 
[T2I Hoi Il6t Il7l IT^ . l24ll . Numerical results for cmissivities from the pair process are taken from Ref. 01 > the photo 
process from Ref. [l!| and the bremsstrahlung process from Ref. [25| . 




FIG. 9: Profiles of temperature versus baryon density ps for varying values of the electron fraction Y e = n e /na, at constant 
entropies per baryon s. 
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FIG. 10: Comparison of neutrino emissivities from the pair, plasma, photo, and bremsstrahhmg processes at constant entropies 
per baryon of s = 4, 6, 8, and 10, respectively. The corresponding temperatures are shown in Fig. 

Traditionally, comparisons of total emissivities from the various processes have been carried out using T and psY e 
as variables. In the context of supernova dynamics, additional insights can be gained by making comparisons along 
constant entropy per baryon adiabats as a function of baryon density ps- We therefore present constant entropy per 
baryon adiabats in Fig. El The results in this figure are based on calculations of the equation of state of matter 
in which the net negative electric charge of electrons and positrons is cancelled by a uniform positively charged 
background of protons, alpha particles, and heavier ions (see, for example, Ref. [26j and references therein). The 
phase structure of matter and the relative abundances of the various constituents including those of dripped neutrons 
at sub-nuclear densities are determined by the minimization of free energy. 

In Fig. I1UI we compare the total neutrino emissivities from the various processes at constant entropy per baryon 
(in units of feg) of s = 4, 6, 8 and 10, respectively. The plasma neutrino cmissivity dominates over the other processes 
considered for entropies and densities that lie at the lower ends of those shown in this figure. Toward the higher ends of 
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entropies and densities, the pair and bremsstrahlung processes are far more efficient than both the photo and plasma 
processes. Once the threshold of T ~ 2m e is crossed, e + e _ annihilation becomes the dominant source of vv pairs, 
the emissivity growing as T 9 . The dominance of the bremsstrahlung process at low entropies and high densities is 
due to the increasing number of participating Pauli-unblocked nucleoli pairs that can benefit from strong interaction 
(pion-exchange) dynamics. Note that the electron fraction Y e has a significant effect only for the pair emissivity, 
which depends directly on the number of e + e~ pairs, which grows with increasing s (or T at a given density). The 
plasma and photo processes are much less sensitive to Y" e , primarily because of the presence of electrons even without 
positrons. The insensitivity of the bremsstrahlung process to Y e is due to the preponderance of nn pairs under these 
physical conditions. 

The relative importance of the plasma process depends on whether or not favorable conditions are encountered 
during the course of supernova dynamics. Simulations peformed to date [53, 13 P om t to the low density (pg < 
10 11 g cm~ 3 ) and low entropy (s < 6) regions encountered a few tens of milliseconds after bounce. For electron 
neutrinos and antineutrinos, the pair reactions are only corrections compared to the beta decay and capture processes. 
For fj, and r neutrinos, densities of ps < 10 11 g cm~ 3 lie below the neutrinosphcrc located at ps ~ 10 12 g cm~ 3 . 
The largest pair rate at the latter density is 5-10 orders larger than the plasma rate in its most relevant region. 
At s = 4, plasma neutrinos are, however, competitive with those from pair and bremsstrahlung processes at the 
neutrinosphcric density (but such low entropies in current post bounce models occur only at much higher densities). 
In order to judge the relative importance for neutrino transport, consult the discussion of thermalization depths for 
the different processes in Ref. [29l |. It must be emphasized, however, that the values of s,T and pg attained in such 
simulations are predicated both by the input microscopic physics and the complex interplay of neutrino transport 
with hydrodynamics (in many cases coupled with convection, rotation and magnetic fields), both of which arc being 
currently improved. The differential cmissivities presented in this work are a part of such improvements in the area 
of microphysical ingredients and their relative importance remains to be ascertained in the future. 

VII. SUMMARY AND OUTLOOK 

In summary, 

• We have calculated the differential rates and cmissivities of neutrino pairs from an equilibrium plasma for baryon 
densities ps < 10 14 g cm -3 and temperatures T < 10 11 K. We have checked that the new analytical expressions 
for the differential cmissivities yield total cmissivities that are consistent with those calculated using independent 
methods. 

• We have developed new analytical expressions for the total cmissivities in various limiting situations. These 
results help us to better understand qualitatively the scaling of the results with physical quantities such as the 
chemical potential and plasma frequency at each temperature and density. A comparison with other competing 
processes, such as e + e~ annihilation, e~7 interaction, and nucleon-nucleon bremsstrahlung, shows that the 
plasma process is the dominant source of neutrino pairs in regions of high degeneracy. 

• Using our results for the differential rates and emissivities, we have calculated the production and absorption 
kernels in the source term of the Boltzmann equation employed in exact, albeit numerical, treatments of multi- 
energy neutrino transport. We have also provided the appropriate Legendre coefficients of these kernels in forms 
suitable for multi-group flux-limited diffusion schemes. 

The significance of neutrino pair production from the plasma process in detailed calculations of core-collapse su- 
pernovae in which neutrino transport is strongly coupled with hydrodynamics remains to be explored. The detailed 
differential information provided in this work may also be of utility in better understanding the evolution of other 
astrophysical systems such as the red giant stages of stellar evolution, the cooling and accretion-induced collapse of 
white dwarfs, burning phenomena atop neutron stars, the explosive stages of type la supcrnovac, and gamma-ray 
bursters. 
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APPENDIX A: EFFECTIVE COUPLING 



The decay of a massive photon or plasmon into a vv pair requires the coupling of the intermediate off-shell e + e~ pair 
to the outgoing neutrinos. Since neutrinos can belong to any of the three flavors, and the loop calculation is performed 
including the e + e~ pair only, different v flavors will have different rates of emission due to different couplings to the 
e + e~ pair. 

The electron-positron couplings to the neutrino-antineutrino pair described by the two one-loop Feynman diagrams 
in Fig. ^are given by the matrix elements: 

M z = I u v — — 7 M (c y - c A j 5 )v u \ (-t) g Q2 _^ /2 — - \v e — - — 7 (cy — c A j 5 )u e j , (Al) 

2 



where in Eq. (|A1|) the superscript i denotes neutrino flavor. The coefficients cy and ca are: 

cy — — and ca = — for v e , and v T 

cy = — — + 2 sin 2 and ca = — — for e, /i, and t. (A3) 

While the exchange of the Z— boson produces all three neutrino and antineutrino flavors (f e>)liT ), W— boson exchange 
produces only v e v e pair. This gives rise to different production rates for different v types so that the electron neutrino 
will have a different emissivity compared to those of the and v T . 

It is economical to rephrase these two diagrams into an effective form dependent on the v type. After performing 
Fierz rearrangement and discarding the Z and W momenta {M% w 3> q 2 ), we obtain the effective e + e~ coupling 

M efJ = ^(m"(1-75K)(^(4-^a75K), (A4) 

f f 

where the couplings C v and C A now depend on the outgoing neutrino flavor. Their numerical values are 

I „ I 

Cv = ^ + 2sm 0w, Ca = r for ^ 

Cy = + 2sin 2 #w, Ca = — ^ for v a and v T . (A5) 

The reduction of the interaction to a single effective interaction significantly simplifies calculations of rates and 
emissivities from the process 7* — > v + v. 
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